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Abstract A direct transient growth analysis for two dimensional, three component perturbations to flow past 
a periodic array of T- 106/300 low pressure turbine fan blades is presented. The methodology is based on a 
singular value decomposition of the flow evolution operator, linearised about a steady or periodic base flow. 
This analysis yields the optimal growth modes. Previous work on global mode stability analysis of this flow 
geometry showed the flow is asymptotically stable, indicating a non-modal explanation of transition may 
be more appropriate. The present work extends previous investigations into the transient growth around a 
steady base flow, to higher Reynolds numbers and periodic base flows. It is found that the notable transient 
growth of the optimal modes suggests a plausible route to transition in comparison to modal growth for this 
configuration. The spatial extent and localisation of the optimal modes is examined and possible physical 
triggering mechanisms are discussed. It is found that for longer times and longer spanwise wavelengths, a 
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separation in the shear layer excites the wake mode. For shorter times and spanwise wavelengths, smaller 
growth associated with excitation of the near wake are observed. 

Keywords Transient growth ■ Stability • Global Modes 

PACS First • Second • More 



1 Introduction 

Hydrodynamic stability is typically studied by the method of linearisation and subsequent modal analysis [ 16]. 
This approach considers the asymptotic behaviour of small perturbations to a steady or time -periodic base 
flow. Such asymptotic behaviour is determined by the eigenvalues of a linear operator arising from the analy- 
sis, describing the time-evolution of the eigenmodes. 

Many canonical problems, such as flow in a channel, permit such stability analysis to be performed about 
a velocity field which depends on a single coordinate. However in more complex geometries we can extend 
the classical hydrodynamic stability analysis to use fully resolved computational stability analysis of the flow 
field IT71 I40I . This is referred to as biglobal stability analysis II35II or direct linear stability analysis in analogy 
to direct numerical simulation (DNS). This approach is able to resolve fully the base flow in two or three 
dimensions and to perform a stability analysis with respect to perturbations in two or three dimensions. This 
methodology does not need to resort to any approximations beyond the initial linearisation and the imposition 
of inflow and outflow conditions. 

In particular the biglobal stability analysis method allows us to consider flows with rapid streamwise vari- 
ation in two spatial dimensions such as the case of interest, flow over a low pressure turbine blade. By postu- 
lating spanwise homogeneous modal instabilities of the form: u'(x, y,z,t) = u(x, y) exp(z'/3z + Xt), asymptotic 
instability analysis becomes a large scale eigenvalue problem for the modal shape u and eigenvalue A. This 
permits use of algorithms and numerical techniques which provide the leading eigenvalues and eigenmodes 
for the resulting large problems, typically through iterative techniques such as the Arnoldi method [40]. This 
approach is extremely effective at determining absolute instabilities in many complex geometry flows, both 
open and closed [7 .,17,15,9,37,36 ,34 22 12,38] including weakly nonlinear stability 1 40 1,231. 

Direct linear stability analysis has not been routinely applied to convective instabilities that commonly 
arise in open domain problems with inflow and outflow conditions. One reason is that such flows are not 
typically dominated by modal behaviour, but rather by significant growth of transients that can arise owing 
to the non-normality of the eigenmodes. A large-scale eigenvalue analysis is not designed to detect such 
behaviour, although for streamwise-periodic flow, it is possible to analyse convective instability through direct 
linear stability analysis QUI . 

To examine this situation, hydrodynamic stability analysis has been extended to cover non-modal stability 
analysis or transient growth analysis [13 39 31 32 ]. This approach poses an initial value problem to find the 
linear growth of infinitesimal perturbations over a prescribed finite time interval. Much of the initial focus in 
this area has been on large linear transient amplification and the relationship of this to subcritical transition to 
turbulence in plane shear flows H20II131 . 

This approach was recently employed in backward-facing-step flows [29,101. With different emphasis 
from the present approach, Ehrenstein & Gallaire [18] have directly computed modes in boundary-layer flow 
to analyze transient growth associated with convective instability and Hcepffner et al 11241 investigated the 
transient growth of boundary layer streaks. A comparable study to the present work [1], studied steady and 
periodic flows past a cylinder. A preliminary study along the lines of the current work, concentrating on a 
steady base flow, is described in 11331 . 

Methods developed for direct linear stability analysis of the Navier-Stokes equations in general geome- 
tries have been previously described in detail [40], and extensively applied [7,6 34,9 11 12,19]. Subse- 
quently, large-scale techniques have been extended to the transient growth problem. In [5] a method suitable 
for such direct optimal growth computations for the linearized Navier-Stokes equations in general geome- 
tries was described in detail. The extension to periodic base flows was presented in [8], for the case of a 
stenotic/constricted pipe flow. The approach has been applied to steady flows past a low pressure turbine 
blades 1331 . a backward facing step HOI and a cylinder, IT), and is the method adopted in this study. 

Recent work [2 3 4] on the flow past the same T- 106/300 low-pressure turbine blade (LPT) as used in 
this study, concentrated on a biglobal stability analysis in order to understand the instability mechanisms 
in this class of flows. At a Reynolds number of 2000, the base flow displays periodic shedding. The work 
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imposed periodic boundary conditions which imply synchronous shedding from all blades. Relevant to the 
current study, the work found that these periodic boundary conditions caused the marginally stable flow to 
go very marginally unstable (Floquet multiplier fi just over 1) at a Reynolds number of 2000 and span wise 
wavelengths L z of approximately L z = 3c to L z = 12c, where c is the projection of the chord length on the 
streamwise axis. This instability is understood to be due to strict enforcement of the periodic boundary condi- 
tions, which is arguably not physical. Relaxing the strict periodicity by using a double-bladed mesh resulted in 
stable eigenmodes (jj, < 1) for all Reynold numbers explored, presumably because small subharmonic effects 
were sufficient to supress synchronous shedding and allow asynchronous shedding. The same geometry is 
analysed here. For our purposes, the flow is best characterised as marginally unstable with Floquet multiplier 
(jj. = 1.00). This work highlights that modal analysis with secondary instabilities does not explain transition 
in this flow. 

Abdessemed et al [2,3 ,' 4j also considered the transient growth problem using a steady base flow at a 
Reynolds number of 895, where significant transient growth up to order 10 5 was observed. The present study 
extends this to periodic base flow at a higher Reynolds number of 2000, where the base flow is periodic. 

The paper is outlined as follows. In Section[2]we outline the direct stability analysis method and the direct 
transient growth analysis needed to determine the peak growth and the associated perturbations. In Section[3] 
we present results for transient growth about a periodic base flow, considering in turn variation by spanwise 
wavelength, period and starting point in the base flow phase. 



2 Methodology 

The flow over the blade is governed by the incompressible Navier-Stokes equations, written in non-dimensional 
form as 



where u(x,f) = [u,v,w](x,y,z,t) is the velocity field, p(x,t) is the kinematic (or modified) pressure field and 
Q is the flow domain illustrated in Figure Q] In what follows we define Reynolds number as Re = U^c/v, 
with Uoo being the inflow velocity magnitude, c the projection of the axial blade chord on the streamwise axis, 
and V the kinematic viscosity. Thus, we non-dimensionalise using c as the length scale, Uco as a velocity scale, 
so the time scale is then c/U^. In the present work all numerical computations of the base flows, whose two 
and three-dimensional energy growth characteristics we are interested in, will exploit the homogeneity in z 
and require only a two-dimensional computational domain. 

We first consider a base flow about which we wish to study the linear stability. The base flows for this 
problem are two-dimensional, time-dependent flows that obey Equations Q] with u = U and P is defined as the 
associated base-flow pressure. The boundary conditions imposed on £2 in the base flow equations are uniform 
velocity U„ at the inflow, fully developed (d„u = p = 0) at the outflow, periodic connectivity at the lower and 
upper boundaries and no-slip conditions at the blade surface. 

Our interest is in the evolution of infinitesimal perturbations u' to the base flows. The linearized Navier- 
Stokes equations governing these perturbations are found by substituting 

u = U + eu' and p = P+ep', (2) 

where p' is the pressure perturbation, into the Navier-Stokes equations and keeping the lowest order (linear) 
terms in £. The resulting equations are 



(u-V)u- Vp + Re~ l V 2 \x in Q, 



(la) 



Vu = inX2, 



(lb) 



(U-V)u' - (u'-V)U - Vp' + Re~ l V 2 u' in Q, 



(3a) 



V u =0 in Q. 



(3b) 



These equations are to be solved subject to appropriate initial conditions and the boundary conditions. 
The initial condition is an arbitrary incompressible flow which we denote by Uo, i.e. u'(x,t — ?o) = Uo(x). 
The boundary conditions we consider are homogeneous Dirichlet on all boundaries, i.e. u'(dQ,t) = 0. As 
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discussed in [10 5], such homogeneous Dirichlet boundary conditions simplify the treatment of the adjoint 
problem because they lead to corresponding homogeneous Dirichlet boundary conditions on the adjoint fields. 

We note that the action of Equations (O (a) and (b) on an initial perturbation u'(x,fo) over time interval T 
may be stated as 

u'(x,T) = ^(T)u , (x,f ). (4) 



2.1 Linear asymptotic stability analysis 

The modal decomposition of this forward evolution operator stf (t) determines the asymptotic stability of the 
base flow U. In this case the solution is proposed to be the sum of eigenmodes, 

u'(x,f) = ^exp(/L/f)u 7 '(x) +c.c, 
and we obtain the eigenvalue problem 

ssl{z)VLj = Ufa, li j = exp(AyT). (5) 

Since for the case of interest U is r -periodic, we set % = T and consider this as a temporal Floquet 
problem, in which case the /J, are Floquet multipliers and the eigenmodes of si (t) are the ^-periodic Floquet 
modes Oy(x,f + T) = Uj(x,t) evaluated at a specific temporal phase. 



2.2 Optimal transient growth/ Singular value decomposition 

Our primary interest is in the energy growth of perturbations over an arbitrary time interval, T. We treat % and 
?o as parameters to be varied in this study. As is conventional 0311 we define transient growth with respect to 
the energy norm of the perturbation flow, derived from the L2 inner product 

2E(u') = (u',u') = [ u'-u'dV, 

JQ 

where E is the kinetic energy per unit mass of a perturbation, integrated over the full domain. The transient 
energy growth over interval T is 

E(r)/E = (u'(t),u'(t)) 

= (^(T)u'(f )X(T)u'(fo)) 
= (u , (f )X*(T)^(T)u'(f )), 

where we introduce srf*{%), the adjoint of the forward evolution operator in (0}. The action of s/*{x) is 
obtained by integrating the adjoint linearized Navier-Stokes equations 

-d t u* = -(U-V)u* + (VU) T -u* - Vp* +Re~W 2 u*, (6a) 



V u* = in Q (6b) 

backwards in time over interval T. The action of the symmetric component operator stf*{T)stf{r) on u' is 
obtained by serial time integration of jz/ (t) and srf*(r), i.e. we first use u'(0) to initialise the integration of 
((3) forwards in time over interval T, then use the outcome to initialise the integration of © backwards in time 
over the same interval. 

The optimal perturbation is the eigenfunction of s^*(x)s^ (t) corresponding to the compound operator's 
dominant eigenvalue, and so we seek the dominant eigenvalues Ay and eigenmodes V; of the problem 

^*{T)^{x)\j=kj\j. 

We use G(t,?o) to denote the maximum energy growth obtainable at time T from initial time to, while the 
global maximum is denoted by G max (fo) = max T G(r,fo)- 
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x/c 

Fig. 1 Hybrid spectral/ftp macro element mesh made up of approximately 2000 elements (left) full mesh, (inset) enlarged view 
around the blade and (right) the double blade configuration. Within each element drawn, a polynomial expansion of order p = 6 
is applied 



Specifically, 

G(Vo) = max(Ay) = max (7) 
./ u'(o) (u'(fo),u'(fo)) 

We note that the eigenfunctions V; correspond to right singular vectors of operator while their 

(L2 -normalised) outcomes Uy under the action of si (t) are the left singular vectors, i.e. 

si(z)\j = OjUj, (8) 

1 /2 

where the sets of vectors Uj and Vj are each orthonormal. The singular values of si (t) are Oj = A. - , where 
both Oj and Ay are real and non-negative. 

While long-time asymptotic growth is determined from the eigenvalue decomposition (0, optimal tran- 
sient growth is described in terms of the singular value decomposition ([H]). Specifically, the optimal initial 
condition and its (normalised) outcome after evolution over time T are respectively the right and left singular 
vectors of the forward operator si (t) corresponding to the largest singular value. The square of that singular 
value is the largest eigenvalue of si* si and is the optimal energy growth G(x,to). 

As already mentioned for an open flow, the most straightforward perturbation velocity boundary condi- 
tions to apply on both the inflow and outflow are homogeneous Dirichlet, i.e. u' = 0, for both the forward 
and adjoint linearized Navier-Stokes equations. The primitive variable, optimal growth formulation adopted 
in this work is discussed in further detail in |[10] |5l and follows almost directly from the treatments given by 
lfl4l l28 241 for strictly parallel or weakly non-parallel basic states. 



2.3 Time integration and spatial discretisation 

Spectral//jp elements 11271 are used for spatial discretisation, coupled with a Fourier decomposition in the 
homogeneous direction. Time integration is carried out using a velocity-correction scheme H26II251 . The same 
discretisation and time integration schemes are used to compute base flows, and the actions of the forward 
and adjoint linearised Navier-Stokes operators. The base flows are pre-computed and stored as data for the 
transient growth analysis in the form of 32 time-slices. The base flow over one period of the evolution is 
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Fig. 2 View of the periodic base flow at (uppermost) the beginning, (second) a quarter of the way through, (third) halfway 
through, (last) three quarters of the way through one base flow period. The figures show spanwise vorticity. 



reconstructed as required using Fourier interpolation. As a check, the results for x/T =4 and ft = were 
repeated with 16 time slices and found to differ by under 0.1%. 

Figure Q] shows the computational domain for the T- 106/300 low pressure turbine blade. The blade ge- 
ometry is approximated by a cubic B-spline interpolation over 200 points to give a smooth flow surface. The 
hybrid mesh consists of approximately 2000 elements, 270 structured elements for the boundary layer around 
the blade surface and an unstructured mesh for the remainder of the field. The elements each have a polyno- 
mial order p of p — 6, with (p + l) 2 degrees of freedom for the quadrilateral elements and (p+ l)(p + 2)/2 
for the triangular elements. Comparison with higher polynomial-order results showed that computations at 
the chosen polynomial order are sufficient to resolve the eigenvalue to about 0.2%. 

To calculate the base flow around which we apply perturbations, a two-dimensional DNS was performed. 
The two-dimensional periodic base flow is calculated at a Reynolds number Re = 2000. The present study 
considers two and three-dimensional linear perturbations to this base flow. 

Spatial periodicity between the planes has been assumed when computing the flow around one single 
blade. The imposition of zero Dirichlet boundary conditions on the outflow for all perturbation computations 
implies a natural limit to the integration period. Care was therefore be taken that the solution does not depend 
on the boundary conditions, in particular the outflow conditions. This was ensured by comparing results using 
an extended domain of twice the length to ensure independence of boundary conditions and solution for the 
cases under study. Doubling the domain length (for the j5 = case) changed the result for x/T = 4 by less 
than 0.2%, for x/T = 8 by less than 2% and for x/T = 16 by less than 8%. All these results are within plotting 
accuracy where presented. In the case of T = 16, the results for the extended domain were presented. 
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Fig. 3 Iso-curves showing transient growth in terms of Fig. 4 Iso-curves showing transient growth in terms of 
/ogio(G) with varying t and L,, with fixed to = 0. /ogio(G) with varying t and /3, with fixed to = 0. 




O 50 



40 



30 
10 




10° 



K) 1 

LJc 



Fig. 5 Maximum growth as a function of L- for t = 0.257\ to = „. , • t , r rrr ^ r, 

r. ^ * j j • i , , ■ ,, . Fig. 6 Maximum growth as a function of L, for t = T, to = 

0. The first and second singular values (squared) are shown. ° * 



3 Behaviour of perturbations to a periodic base flow 

Throughout the study the Reynolds number is fixed at 2000, corresponding to a periodic base flow (after the 
first Hopf bifurcation at Re=905). The shedding period of the base flow is T = 0.2676 and the phase of the 
initial condition in the base flow shedding cycle is fixed at to = as depicted in Figure [2] and as defined in the 
previous section, unless otherwise specified. As indicated in the introduction, the flow is best characterised as 
marginally unstable with Floquet multiplier (fi = 1.00). 

The dependence of growth on z and L z is shown in Figure [3] and the same data is shown for j3 = 2n/L z 
in Figure [4] We find that, while for spanwise wavelengths above about 10c, the growth achievable depends 
little on spanwise wavelength L z , below L z ~ 10c, the optimal growth in a given T is limited, particularly for 
longer times. As L z tending to infinity represents a purely two-dimensional problem, it appears that it is the 
two-dimensional perturbations that dominate potential instabilities at this Reynolds number. Growth of about 
over 10 4 is achievable for T ~ 107\ For longer times the disturbance has convected too far downstream to be 
of interest. 

To examine the case of shorter T optimal disturbances we take z = 0.25T and vary L z as a parameter. 
Figure U] shows the first two leading growth values associated to the two most significant optimum modes. It 
can be seen that the growth is moderate and concentrated at short L z . The maximum occurs at Lz = 0. 125c for 
this integration time. For slightly longer times (t = T and above) the growth rises with L z and then flattens 
out, as shown in figure [6] This effect is also evident in figure [7] 

Taking the peak of the first mode at T = T (see Figure [6} and varying T produces the results shown in 
Figure [8] Again, strong growth is seen up to T ~ AT after which a plateauing is observed as the disturbance 
aligns to the least stable eigenmode (nearly marginally stable). The optimal mode for z = T and L z = 1.95c 
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Fig. 8 Maximum growth as a function of T for L, = 1.95c, 
Fig. 7 Maximum growth as a function of L z for % = 16T, to = t Q = 




Fig. 9 Normalised leading mode initial condition (left upper) and final state (left lower) at L z = 1.95c, for z = T, at zero initial 
phase, showing spanwise vorticity. The right-hand plot shows the same for the mode associated with the second singular value. 



is shown in Figure [9] which shows a disturbance beginning at the trailing edge and the shear layer exciting 
the near wake. The spanwise constant optimal mode for T = 87/ is shown in Figure[TT] In contrast to the short 
T mode shown in the left panel of figure [9] the wake mode is excited some distance downstream. Although 
the initial disturbance also involves separation at the shear layer and a disturbance at the trailing edge, the 
initial disturbance extends further up the blade surface. Figure [lOl shows iso-surfaces of streamwise vorticity 
at (O x = 0.6 of the same mode combined with 10 x the base flow. 

The preceding discussion has assumed an initial phase of fo = (relating to the point in the shedding 
cycle) at t = as defined in figure [2] An investigation was also carried out to find the effects on the maximum 
attainable growth of varying this initial phase. At short growth horizons T, the growth achieved is only weakly 
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Fig. 11 Normalised leading mode initial condition (upper) and final state (lower) at /3 = (2D case), for % = 87\ at zero initial 
phase, showing spanwise vorticity 
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Fig. 12 fie/y G T with varying to and /3, t = 0.25r; ("n'gtoj /ogio(G(r)) with varying % and ?o, /3 = 



dependent on L z and is achieved just before ?o ~ 0.2 as can be seen in figure [12] Furthermore the figure also 
shows the dependence on initial phase remains weak over a range of integration times T. We conclude from 
the figure that the dependence of growth on the initial phase ?o is relatively unimportant. 

The relevant previous study [4 | demonstrated that although subharmonic effects were small, they were 
sufficient to induce asynchronous shedding and prevented the flow from going unstable. For the present prob- 
lem we do not consider the asymptotic or long-time behaviour of perturbations but their behaviour over a 
shorter time horizon. 

Results for the double-bladed mesh were tested for times up to eight shedding periods and agree well 
with the single blade situation. At these longer times, consideration of the least stable eigenmode would be 
more suitable. This leads to the conclusion that subharmonic effects are relatively unimportant in the transient 
growth problem and the assumed periodic boundary condition is a valid means to reduce the computational 
domain for the problem under study. 



4 Discussion and conclusion 

The transient behaviour of perturbations to linearised flow past a periodic array of T- 106/300 low pressure 
turbine fan blade was investigated. The analysis was carried out at a Reynolds number of 2000 associated 
with periodic vortex shedding, used as a periodic base flow. 

It is known from asymptotic analysis [4 1, the flow past an array of these turbine blades is marginally stable. 
The current analysis shows that long wavelength optimal perturbations associated with long time integration 
periods convect far downstream and eventually align with the asymtotically least stable eigenmode. The dis- 
covery of converging optimum growth for long integration times confirms the lack of a strong asymptotic 
instability. 

It is found that the long-wavelength perturbations tend toward a purely two-dimensional case and that 
these perturbations are associated to maximum optimum growth in the asymptotic case approached by long 
time-integration. However, as we already know from two-dimensional DNS, the two-dimensional baseflow 
is stable for nonlinear flows and therefore also stable in a linear sense. It therefore may be assumed that the 
identified optimum growth associated to long wavelength perturbations is less significant than perturbations 
that are limited in spanwise wavelength. When considering short integration times it is indeed found that 
short-wavelength perturbations are of higher significance. 

Short integration times have maximum optimal growth at shorter spanwise wavelengths and convect only a 
short distance downstream, exciting the near wake. We might hypothesise that in the presence of the neglected 
nonlinearity, the associated optimum modes cause shear layer separation, triggering wake instability. This 
would have to be demonstrated in a full nonlinear DNS. Furthermore, the spanwise length of a real LPT blade 
is naturally limited, giving an additional reason why short wavelength perturbation are more important than 
long wavelength disturbances of theoretical spanwise extent. 
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Our results are consistent with the understanding that transient growth mechanisms are associated with 
shear in the base flow feeding the perturbation energy growth. The results presented are consistent there- 
fore with the cylinder results reported in the literature (both on transient growth mechanisms lfl"l l331 and on 
receptivity via the adjoint [21]) and with the previous work on the LPT fan blade used in this study 

It is hoped that the understanding developed will prove useful in controlling laminar boundary layer sepa- 
ration with a view to improving performance. For instance, the spatially periodic pattern in the shear layer can 
be quantified in the time-domain indicating disturbance frequencies susceptible to optimum amplification. 
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